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Let H^'^^ = (g) , . . (g, }j\<i-n] tensor product of Hilbert spaces and let ro be the closest 

separable state in the Hilbert-Schmidt norm to an entangled state po- Let fo denote the closest 
separable state to po along the line segment from I /N to po where I is the identity matrix. Following 
|lj a witness Wo detecting the entanglement of po can be constructed in terms of J,to and fo. If 
representations of to and fo as convex combinations of separable projections are known, then the 
entanglement of po can be detected by local measurements. Giihne et. al. in Q obtain the minimum 
number of measurement settings required for a class of two qubit states. We use our geometric 
approach to generalize their result to the corresponding two qudit case when d is prime and obtain 
the minimum number of measurement settings. In those particular bipartite cases, ro = fo. We 
illustrate our general approach with a two parameter family of three qubit bound entangled states 
for which ro 7^ fo and we show our approach works for n qubits. 

In ^ we elaborated on the role of a "far face" of the separable states relative to a bound entangled 
state Po constructed from an orthogonal unextendible product base. In this paper the geometric 
approach leads to an entanglement witness expressible in terms of a constant times / and a separable 
density po on the far face from po. Up to a normalization this coincides with the witness obtained 
in [0] for the particular example analyzed there. 



I. MOTIVATION AND NOTATION 



An important question for quantum information theory is how to determine if a given state is entangled. Physically, 
one would like to do this using local measurements and classical communications. Testing for entanglement is closely 
related to Bell inequalities J || and subsequent elaborations of Bell's inequalities Recently other tests have been 
suggested, such as that in p ] which relies on the theory of positive operators and on eigenvalue estimation. 

An alternate approach, which is experimentally realizable, is to define local correlated measurements motivated by 
some knowledge of the structure of p itself, and this approach has been elaborated in To describe the problem, 
we first define the mathematical context. Specifically, we assume we are working with n distinct systems so that p 
is represented as an A^ x A^ density operating on the tensor product Hilbert space = ijl'^il (g) • • • (g) ffl''"]. The 

set D of such N x N densities operating on is a compact convex subset of the real Hilbert space M of N x N 
Hermitian matrices where the inner product is defined by {A, B) ^ Tr [A^B] . (S ince the matrices are assumed to 
be Hermitian, the notation "t" denoting the Hermitian conjugate appears to be redundant. However, we will have 
occasion to use the inner product for more general matrices.) The set of separable densities S is defined as the convex 
hull of the separable projections tti (g • • • (g 7r„ , where tt^ is a projection on H^'^''\ Since is a compact convex subset 
of D one can test for entanglement by showing p is separated from S* by a hyperplane in M . Geometrically the idea 
is clear. Mathematically it reduces to finding a Hermitian matrix W with the property that Tr [W p) < Q < Tr (Wa) 
for every density a in S. The existence of such a is guaranteed by the general theory of convex sets in Hilbert 
spaces, and W is known in the quantum information literature as an "entanglement witness" . A nice introduction to 
the subject and an overview of some of the literature can be found in |^. 

In the context of two qubits Giihne et. al. in ||^ assume the general form of a two parameter family of densities 
p which includes a maximally entangled state pq. They construct an entanglement witness using the eigenvector of 
the partial transpose of p with the minimal (negative) eigenvalue and find that the resulting witness does not depend 
on either of the parameters. Since the separating hyperplane contains a face of the separable states, it is optimal in 
the sense that no witness detects a strictly larger set of entangled states. (See |^ for the definitions and 0] for an 
exposition related to the approach used in this paper.) 

In the authors showed how an entanglement witness Wq sensing an inseparable po can be constructed if one 
also knows the nearest separable state tq: 

IIpo - Toll = inf {||po - (t\\ : (J e S} . 

Since the norm is a continuous function and the set of separable densities is compact, tq exists, although actually 
computing it is not an easy problem in general. The entanglement witness is defined by 
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Wq = To + CqI - Pq, (1) 



where I is the N x N identity matrix and 

Co = Tr{To (po - To)) . 

Details and examples of this construction are given in where it is shown that Wq is linked to the geometry via the 
induced inner product 

{{po - To) , (p - To)) = Tr {{po - To) {p - To)) = ~Tr [Wop] . (2) 

In particular the separating hypcrplane contains the "nearest" face of S consisting of separable states a such that 
a — To is orthogonal to po — To- Equation (H) can be used to show that the extreme separable projections in the convex 
representation of tq must lie in the hyperplane, and that if any separable a in the nearest face has full rank then Wo 
is optimal. 

It was first shown in [ pO[ that there is a neighborhood of the normalized identity or completely random state, 
Do = {1/N)I, in which every state is separable . Given that fact, it follows from another compactness argument that 
there is a nearest separable density to po along the line segment [Do, po]'- 

To = {I - So) Do + sopo (3) 

with < So < 1. While Tq and tq differ in general, in certain examples they are the same which simplifies the analysis. 
Thus we have the following general result. 

Theorem 1 Suppose po is inseparable. Using the notation above, the Hermitian matrix 

Wo = I ( Co + ^ ) + To - —fo 
V Nso J so 

is an entanglement witness for po and is optimal if the nearest face contains a separable density of full rank. □ 

Thus if one knew the convex representations of To and fo in terms of tensor products of local projections, one 
could define specific coordinated local measurements that would experimentally detect the entanglement of po via 
Tr (poWo). Finding tq and fo is in general difficult but can be done in a variety of special cases. The examples we 
present include those analyzed in |^ as well as a two parameter family of three qubit bound entangled densities for 
which To and fo differ. (A bound entangled state is entangled but has positive partial transposes.) 

Another result in |2j is that three sets of coordinated local measurements is the minimum number required in their 
two qubit context and an explicit representation of the three measurements was given. It was also asserted that at 
least d + 1 such measurements would be required for a corresponding d x d system, but no suggestion for achieving 
that bound was provided. We show how the geometric approach to entanglement witnesses provides a unifying theme 
and leads to a concrete construction for the d x d case when d is prime. 



II. TWO QUBITS 

In the two qubit case, the use of the nearest separable density clarifies some of the methodology and suggests the 
generalization to the d x d case. Following ref. we take 

P = PPa + (1 - p) <J. (4) 

Pa is defined by the state a |00) +b where a and b are real with + 6^ = 1, p is a parameter between and 1, and 
cr is a density close to the normalized identity, \\a — Do\\ < 6. The density a represents noise that is close to Do, the 
completely random state. The idea is to define a separating hyperplane Wo based on po = | (|00) + |11)) ((00| + (11|) 
and investigate what inseparable states p are detected by Wq. 

It has been shown in a number of places that the closest separable state to po is 

2 1 

To = -^Do + -po, (5) 

so that the roles of tq and fo coincide. (References and details are given in [Q.) One computes cq — Tr (tq {po — To)) = 
^, and then eq. (0) gives 
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/ 0-1 

^° = 3 1 1 ■ 
V -1 

This differs from the optimal witness found in Q only because of the use of a different Bell state and is a special case 
of the general theorem above. 

As an application we have the following result. 

Lemma 1 A sufficient condition that p — ppa + (1 — p)cr, where \\a — Do\\ < S, is not separable is that 

1 + 4(5 



Aab+1 + 4(5 



< P- 



Proof : 



Tr (Wop) - pTr (WoPa) + ((1 - p)) Tr (WoDo) + {I ~ p) Tr (Wq (a - Dq)) 
+ ^{l-p) + il-p)Tr{Wo{a~Do)) 



3 6 
-2abp 1 , , ,2(5 

where we have used the Cauchy-Schwarz inequality in the last step. Setting the final expression to be less than 0, we 
obtain the desired inequality. Note that if 6 — and a = l/\/2 we obtain the well known sufficient condition 1/3 < p 
for inseparability of ppo + (1 ~ p) D^. □ 

Having defined Wq we need to show that the measurement can be effected by three types of coordinated local 
measurements. We combine eq.(|l|) and eq.(^ to obtain 

Wo = h- 2to (7) 

and then use the representation of tq as a convex combination of six separable extreme points in the face of the states 
of S in the separating hyperplane: 



f , I CTq + CTz ^ CTq + CFz\ , / CTO — (To — Cr 



6'V 2 2 y V 2 2 

+ CTa; (To + da; \ /cTq — (7^ _ (Tq — (Ta 



(8) 



Co + CTy CTo — CTjy \ / CTq — CTj, CTo 



cr,, 



Thus one takes coordinated local measurements along the x, y, and z axes of the Bloch sphere to compute Tr (Wop) = 
I — 2Tr (tqp). As shown in j^, this is the minimal number of coordinated local measurements which are required. 

III. THE dxd CASE. 

The approach used above immediately generalizes to the bipartite dxd case when d is prime (d ^ 2): we take an 
entangled "base" state po for which we can compute the nearest separable state tq and thus Wq. We again consider 
the family of densities p — ppa + (1 p) a, where a is close to the state Do, and define 



Pa = \tpa) {lpa\ wkcre |V'a) = X! l^^'* 



k=0 



with real ai~ such that a| = 1. po is the state with Ok = l/Vd, and 
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is the closest separable state. (See for the general result and references.) Again, tq coincides with tq, simplifying 
the problem |p],|l2|-p^. From eq.(|l|) the optimal witness for po is Wq = jj^I — dr^, where Cq — ^j^jqiiy- The problem 
now reduces to finding analogues of the Pauli matrices which can be used to represent tq as an appropriate convex 
combination of projections, as in eq.( ||). Fortunately that analysis already has been done. 

In ref. the authors observed that the (real) Pauli matrices can be viewed as discrete Fourier transforms of four 
"computational" basis matrices. Using an analogous basis for d x d matrices and the corresponding discrete Fourier 
transform, one is able to define (P orthogonal unitary matrices 

Ud = {Su:u=ij,k), 0<j,fc<d} 

where — 5'cn.n} is the d x d identity. (These same matrices had been derived independently and in a different 
manner by Fivel [ |l6[ who used them in a study of Hamiltonians on a discrete state space. He also derived several of 
the properties we include below.) As with the two qubit case, one can define sets of tensor products of projections, and 
it turns out that tq can be written as a convex combination oi d+1 such sets in strict analogy with the representation 
in eq.(^). These d+1 sets of projections correspond to the coordinated local measurements required in ||^ for local 
detection of entanglement of d x d states. 

We briefly summarize the necessary properties of these d-level "spin" matrices and relegate proofs to the Appendix. 
By definition 

d-l 



r=0 



where addition is modulo d and rj = exp (27ri/rf) 

Tr 



'5'o-i,fci)'^02,fc2) 



dSji,j2^ki,k2 



expresses orthogonality, and thus Ud is a basis for d x d matrices. 

Using tensor products of the Ud spin matrices, we find that for prime d 



1 

•^0 = -1— 



1 



k=0 




k(d~k),no + 7. M« 7. ^(kj){kd-kj),kk 



(10) 



where Sij^ — Si^k (X" Sj^i, and it remains to show that each of the ^-summations can be written as a sum of tensor 
products of complete sets of projections. When d is odd and u = {j, k) ^ (0, 0), 

m=0 

is a (Hermitian) projection, and {Pu (r) : < r < d} is a complete set of orthogonal projections. If Uj — {j, 1) and 
Vj = (d — J, 1) for < j < d, then (surpressing the subscript on Uj and Vj) 

^ d-l ^ d-l 

X! ^{k3){kd-kj),kk = j^^Pu{r)<E) Pv (d-r). 

k=0 r=0 

The first summation ( p^ ) has an analogous representation if u = (1, 0) and v — (d — 1, 0). 

This completes the proof: the entanglement witness Wq can be realized in terms of the identity and a separable 
density which in turn can be written as convex combination of d + 1 sums of tensor products of complete (local) 
projections. This attains the lower bound for the number of coordinated local measurements as asserted in 

As in the two qubit case, the entanglement witness Wq detects entanglement for a range of densities of the form 
p = ppa + (1 — p) cr. The computation is similar to that for Lemma 1, and we omit the details. 

Lemma 2 p ~ ppa + (1 — p) cr is inseparable provided 

where \\a — Dq\\ < 5. When d — 2, this reduces to the inequality in Lemma 1. □ 
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IV. A THREE QUBIT EXAMPLE 



The geometric appproach also works for a particular two parameter family of three qubits which have positive 
partial transforms but are inseparable. Since these densities are not generated by complete UPB sets, it is not clear 
that other techniques can be used to define an appropriate entanglement witness. 

Let — 1/8 < c, d < 1/8 and define the three qubit density matrix 



p{c,d) 



/1/8 



1/ 








1/? 

Vl/8 



1/8 



1/8 
d 



l/8\ 



1/ 






c 

d 

1/8 

1/8 

1/ 

1 



It is convenient to identify p{c,d) with the four vector {d,c, 1/8, 1/8) defined by the negative diagonal. We will use 
this notation for densities with analogous structure. Further, it simplifies calculations to use m = (c + d) /2 and 
t = {c — d) /2, and we abuse notation by writing p (m, t) for the same density and (m — t^m + t, 1/8, 1/8) for its four 
vector, 
with equal 



The following result is proved in for analogous densities p{c,d) for n qubits defined by the 2"^^ vector 
ual numbers of c's and d's and 2"-^ntries of 1/2", (d, • • • , d, c, • • • , c, 1/2", • • • , 1/2"). 



Proposition 1 p{c,d) has positive partial transposes and is completely separable if and only if c = d (t = 0). □ 

In the d X d cases analyzed above, the line segment from Dq to pq was orthogonal to the nearest separable face, 
and that property characterized tq, the nearest separable density to po. Unfortunately, as shown in 0], that perpen- 
dicularity is lost when one goes to three systems and the nearest separable state tq to po on the line segment [Dq, po] 
does not coincide with the closest separable state tq. However, one can still take advantage of the geometry provided 
To lies in the nearest separable face to po. 

To pursue this idea for p (c, d), we need some additional notation. Without loss of generality we take c > d so that 
t > 0, and let cri = and (T2 — (^y Let ctq denote the 2x2 identity and define 



jki 



[(To (8) (To <Xl (To ± (Jj 



> (Tfe <Xl (T/ 



(11) 



where j, fc, and I will take the values 1 or 2. It is an easy exercise to represent such a as an average of four 
projections and to confirm that 



p{m,t) 



4m I P+i 



4^ (^212 + ^122, 



1 



2 

- 8tDo. 



4m R 



221 



As it happens, a study of the m = case is key to the analysis, and we take as a candidate for tq (0, t) the 
normalization of the first part of p(0,t): 



fo(0,t) = 



1 



l + 8t 



l{Pl\l+P22l)+^t {P^^^+P{ 



+ \ 
22) 



(12) 



fo(0,i) is obviously separable but not so obviously the last separable state on [I?o, p(0, t)]. We confirm that property 
later. To see if tq (0, i) lies in the same face as fg (0,i), we take normalized combinations of the four Pj^j in the 
equations above and minimize the distance to p (0, i) , finding the separable density with four vector 



t t I t 1 t 

'2' 2' 8 ~ 2' 8 ^ 2 



Using this for tq (0, t) we find co (0, t) = j and Wq (0, t) = | Wq where 
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These heuristics work splendidly, and we also find that Wq detects the entanglement of all p (c, d) with d < c. This 
is illustrated in fig.l where the separating plane is shown. Of course, the geometry is more complicated because the 
hyperplane is not two dimensional. As t decreases to 0, p(0,t) moves to the center of the line segment [-Pni, -P22i]- 
As t becomes negative (d > c), tq and tq move onto a new plane where -P212 ^^'^ ^122 replace P212 ^122- Recall 
that Do lies at the center of [Pj^f.> ^ijk\- m 7^ is easily visualized. 

Proposition 2 Let d < c. Then tq (to, t) with four vector (m is the closest separable density 

to p{m,t). Wq is an entanglement witness for every density in {p{c,d) : —1/8 < d < c < 1/8}. 

Proof: Set m = and let tt = tti ® 7r2 (8) tts denote any separable projection. Since Tr {Wop {0,t)) < by 
construction, it suffices to confirm that Tr {Wop.) > 0. Defining tt/j — |V'fe)('0fc| where 



COS 6k 
sin Ok 



we obtain 



where 



2Tr (Wop) - 2 - sin (20i) sin (2^2) sin (2^3) C {ip, (^3) 



(13) 



C (Vi 7 <yf2 : '^3) = COS (</?! +ip^+(p3)+ COS {ip^ + - ipa) 

+ COS {ip^ fs) - COS {ip^ ~p^~ Vs) ■ 

The phase angles p^ can take any value while < 6'fc < 7r/2. Confirming that the right side of eq.(p^) is non- negative 
is a familiar Bell-inequality computation and proves the assertion when m — Q. It follows from comments after eq.(|^) 
that To (0, i) has to lie in the separating plane and is thus the closest separable state to p (0, i) along [Z?Oi P (0, t)], 
justifying the notation and the assumption made earlier. 

The generalization to non-zero m is straight-forward. Using the asserted form for tq {m,t), it's easy to check that 
To (to, t) is separable, that p (to, t) — tq (to, t) = p (0, t) — tq (0, t) and also that cq (to, t) — ^. It follows that tq (to, t) 
is the closest separable state to p (m, t) and Wq (to, t) — Wq (0, t), completing the proof. □ 

From p (to, t) = (1 -|- 8t)fo (to, t) — 8^1)0 and the form of tq (m, t) we can express the entanglement witness in terms 
of the identity and explicit separable states: 



Wo (0, t) = t 



2 {P^ii + P22I + -^212 



P+ ) 
^122) 



Again we have shown that local detection of entanglement can be defined using the explicit representations of tq and 
To as convex combinations of separable projections. 



V. GENERALIZATION TO TV QUBITS 



In ref. [|l7| tq was computed for the po generated from the n-qubit GHZ state defined by 

1 



Pq = IV'o)(V'o| where \ipo 
and j = {j , ■ ■ ■ , j) ■ It was shown that 



V2 



(|o> + |i» 
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fo = (1 - sq) Do + sopo (14) 
= soA + (l-so)Q (15) 

where sq = 1/ (2"-i + l), A = i (|6)(0| + and Q is a matrix 2" x 2" matrix with entries 1/2" on the diagonal 

and in the upper and lower corners. It is clear that A is a convex combination of two separable states. In Q was 
expressed in terms of 2"^^ separable states. In ||l| we also computed tq and it can be shown that tq can be expressed 
as a convex combination of A and Q. This is another example of a case when tq ^ tq but both densities lie on the 
near face. Applying Theorem 1 the optimal entanglement witness can be written as 

Wo = aI -bA- cQ (16) 

with a, b, and c positive. In the two qubit case this result reduces to eq.(^. 

VI. FAR FACE CONSTRUCTIONS 

There are cases when an entanglement witness can be defined in terms of the identity and a separable state without 
computing the nearest separable density explicitly. In a technique is described for the construction of inseparable 
densities with positive partial transposes, using orthogonal unextendible product bases (UPB). This clever approach 
assumes a set of m separable orthonormal states B = {\ipk) , ^ k < m} where each \ipk) is a tensor product of 
states in their respective Hilbert spaces and where the orthogonal space B-^ contains no separable projections. If 
^J'k = Ifk) {^k\ and one defines 



m ^ — ^ 

then 



Mo 

TO 

fe=l 



N m 
Po = -T-r Da - — ^0 

I\ — m 1\ — m 

can be shown to be an inseparable density with positive partial transform. A number of examples of orthogonal UPBs 
are given in p8[ | and in subsequent papers such as [ p9[ and [ p O | . The ideas in ||^ also apply in this context and are 
illustrated there using the two qutrit example "TILES" of |18[. 

In some consequences of the geometric structure implicit in this approach are developed. For example, it is 
clear from the equation above that Dq lies on the line segment [/ig, po]- If one denotes by Fq the face of the separable 
densities S containing /io, then, in the context of the real Hilbert space M, Fq is orthogonal to that line. It is shown 
in j|] by a compactness argument that there is a positive e such that 

< -i = inf{Tr[^of^] , <J S} 
m 

and thus that the face Gq = {a ^ S : Tr [fiQcr] = ^} is non-empty. In this context it is shown in ||^ that 

eN 

< so = 1 < 1. 

m 

A consequence of this approach is that a separating witness Wq for po can be defined using (§) with 

To = (1 - So) -Do + soPo- 

In this construction tq is not necessarily separable but is defined by the intersection of [/io,Po] and a hyperplane 
containing the "near face" Gq. Since both po and To can be written explicitly in terms of /^o, which is separable, then 
once e is known we can again express the entanglement witness in terms of the identity / and a separable density 
whose convex representation is known: 

eN / e 



Wo = (mo - ^/) 

iV — TO V m / 



Thus, the required coordinated local measurements are defined explicitly by the original set B and there will be no 
more than to different settings. Geometrically Wq is expressed in terms of the identity and /io, which lives in the far 
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face Fq, on the "other side" of Dq from po. In the special case discussed in 0, this is the same witness as derived 
there, up to a muhiphcative constant. 

Consider separable densities p,b = '^i^PkfJ'k in the face Fq that are also close to /zq. Let b denote the reciprocal of 
the largest of the coefhcients pk ■ Then it is easy to define inseparable densities 

NDo - b^ib 

with positive partial transposes that are on the boundary of the set of densities D and are close to po- Moreover Wq 
can also serve as an entanglement witness for these densities. In fact, using the same notation as above, one can get 
a "frustram" of states of the form 



P= {1 - p) (J + PPb 
which lie in D on the po side of the hyperplane defined by Wq, provided 



pjm-b) 1-p jl ~p)S e 
N-b +^ + ^7^<m 



where ||(T — -Doll < We omit the details, repeating instead that the Euclidean geometry of M provides an extremely 
useful context for examining questions of this sort and that the use of eq.(0) gives a unifying geometric approach for 
constructing entanglement witnesses. 

We should note that the effects of dropping the hypothesis that the states in B are orthogonal is also discussed in 
1^ , and weaker conditions on the states in B are given which allow the construction above of inseparable states to be 
generalized. In particular, one can perturb the orthogonal UPB case, losing orthogonality but preserving enough of 
the structure to allow the analysis to go through. The cost of this generalization, however, is that the resulting states 
do not automatically have positive partial transposes. 



VII. SUMMARY 



In this paper we have used a geometric definition of an entanglement witness Wq detecting an inseparable state po 
to show that Wq always has a representation leading to entanglement detection using coordinated local measurements. 
This approach gives essentially the same witnesses and the same coordinated local measurements as derived in |^ for 
their particular two qubit case. When coupled with the generalized "spin" matrices defined in ||l^, it also achieves 
the lower bound asserted in |^] for the number of coordinated local measurements for the analogous d x d case, at 
least when d is prime. We also illustrated the use of the geometry by applying the methodology to a two parameter 
family of three qubit bound entangled states for which tq and tq differ. The strength of the geometrical approach is 
further illustrated by applying it to the n qubit case. In the case of inseparable densities constructed using orthogonal 
UPBs, the geometric approach also applies, but produces a representation using a "far face" separable density. 

VIII. APPENDIX 

By definition the <S'(j_fc) "spin" matrix is defined as 

d-l 

r=Q 

where addition is modulo d and 77 = exp(27ri/d). If u denotes {j,k), then Su has trace unless u equals e = (0,0). 
Orthogonality follows from 



Tr 



S, 



t 



(il,fcl)'^(i2,fe2) 



Tr 



= Tr 



Similarly, one can calculate some useful relations such as So^iSi,o = r]Si,oSo^i, Sj^k = (5'i,o)'' ('S'0,1) , {Sj^kY^ = 
^jkm{m-i)/2 g^^^^^^ and S^- = rj^^ Sd-j,d-k- Unlike the Pauli matrices, the Su are not necessarily Hermitian, but 
they are unitary and can play a role analogous to that played by the Pauli matrices. 

Any d X d density a can thus be written as a linear combination of these spin matrices, and we have 



where we use u — (j,k) in 



To represent a density such as tq defined in eq.(^) on the tensor product space ijl'^' (8) H^'^\ we use the set of tensor 
products of the spin matrices as an orthogonal basis. A direct calculation or an invocation of eq.(16) of ref. ||l^ gives 



To 



1 

d2 



d-1 d-1 



,kk 



k=0 1=0 



where Sij^ki — St^k ® Sj.i. It is at this point that we require d be prime. Then for given i and k ^ Q there is a unique 
j such that i = jk (mod d), and we can rewrite tq in the form 



1 



To 



d+1 
1 

d+1 



d-1 d-1 



-^SoOfiO ^iid-i),00 +J2^Y1 ^i''- 

1=0 fc=l j=0 

^ d-1 d-1 / ^ d-1 

I2 X! ^k{d-k)m +^\-n^ ^(k. 



(kj)(kd-kj),kk 



d^ 



k=0 



j=0 



^2 / ^ ^(kj){kd-ki),kk 
k=0 



It remains to show that each of the expressions involving a /c-summation is a summation of tensor products of 
projections from a complete set of orthogonal projections. That is, each summation corresponds to correlated local 
measurements, and tq is realized by d + 1 such summations. 

We begin by defining a complete set of projections in terms of the spin matrices, a construction which corresponds 
to that in the spin 1/2 context. 



Lemma 3 (Reference Jjq/j Let d > 2 be prime and let e ^ u = [j, k). Then if 

d-1 d-1 
m—O m—0 

{Pu {r) : < r < d} is a complete set of trace one, orthogonal (Hermitian) projections. 

Proof: Pu (r) has trace one since the only term with non-zero trace is the m = term. From the definition 

d-1 d-1 

Pu (r) Puis) =-J^J2J2 



m—O 71—0 
d-1 d-1 



_ 1 \^ \^ „(m+n)r jklm{m-l)+n{n-l)]/2 jkmn q n{ 
— ~n 2-^ n n n '-'(rn+n)uV 



s — r) 



m—O n—0 

Make the substitution t = m + n in the last expression and collect terms to obtain 

d-1 



Pu (r) Puis)=l'-:Y: V"-'-''^ ^ E V'W^'-'^^'S,U 



\ n=0 / t=0 

thereby obtaining both the orthogonality and Pu (r) Pu (r) — Pu {r) ■ Finally 
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m—O m—0 



d-1 d-1 

\ ^ —inr+rajk+jkm{m—l)/2Q ^""^ nr+j/i;n(n— 1)/2 q 

— ~; / '/ 'J—mu — ~; / '/ '-'rati 

m— U n—i) 

= Pu {r) . 

These steps actually introduce a factor of the form r]^'"^('^~'^y^ which equals 1 for odd integers. However, if d is even 
and j and k are odd, 7^J'='^('*-i)/2 and the proof must be modified. □ 

Having defined complete sets of projections, we are ready for the final technical result. 

Proposition 3 Let Uj — {j, 1) and Vj — {d — j, l)/or < j < d. Then 

^d-i / ^ d-1 \ 

3 E ® {d-r)= H2 E S(k3)(kd-kj),kk ■ 

If X ~ (1,0) and y = (d — 1, 0) , then 

^ d-1 ^ d-1 

Px (r) (E) Py (d - r) ^ —Y^ Sk(d-k}, 00- 

r=0 k=0 

Proof: The proof is just a matter of navigating the notation. Suppressing the subscript, 

d-1 



i ^ P„ (r) ^P^{d-r)^^J2 (^«)' ® (^-)" E '^^'""^'^ 

r=0 fc,n r 



]j E ^ ^^^^b-fe(fc-l)/2+(<i-j)fc(fe-l)/2] 



1 1 '''^ 



^2^"--"- ^2 



'(fej")(fc(i-fcj),fcfe 

as required. The proof of the remaining assertion is similar, and we omit the details. □ 
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FIG. 1. The plane containing ro,fo, and the four separable densities of eq. is illustrated. The separable states lie below 
the plane and the random state -Do is shown. 
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